Admittance of the SU(2) and SU(4) Anderson quantum RC circuits 
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We study the Anderson model as a description of the quantum RC circuit for spin- 1/2 electrons 
and a single level connected to a single lead. Our analysis relies on the Fermi liquid nature of the 
ground state which fixes the form of the low energy effective model. The constants of this effective 
model are extracted from a numerical solution of the Bethe ansatz equations for the Anderson model. 
They allow us to compute the charge relaxation resistance Rq in different parameter regimes. In 
the Kondo region, the peak in Rq as a function of the magnetic field is recovered and proven to be 
in quantitative agreement with previous numerical renormalization group results. In the valence- 
fluctuation region, the peak in Rq is shown to persist, with a maximum value of h/2e^ , and an 
analytical expression is obtained using perturbation theory. We extend our analysis to the SU(4) 
Anderson model where we also derive the existence of a giant peak in the charge relaxation resistance. 

PACS numbers: 71.10.At, 73.63.Ay, 72.15.Qm 



I. INTRODUCTION 

High frequency transport experiments aim to control 
and probe the coherent motion of electrons in real time 
[I}]3]. Continuous technological progress have paved the 
way for the study of dynamical properties of mesoscopic 
systems. A typical example is the quantum RC circuit 
[5H7], where a quantum dot is connected to the quan- 
tum Hall edge states of a two-dimensional electron gas 
(2DEG) through a quantum point contact (QPC). As il- 
lustrated in Fig. [T] the quantum dot forms a mesoscopic 
capacitor with a top metallic gate and the charge in the 
dot can be changed periodically by applying an AC drive 
through the gate voltage Vg. For large gate voltage mod- 
ulations, this device acts as a single electron emitter |S1- 
[T^ and the effect of sweeping the last occupied level of the 
dot across the Fermi energy has received a broad theo- 
retical attention [13 16 . The quantum RC circuit is also 
promising for efficient charge readout in a quantum dot 
device |17j or to detect topological excitations [T51 [TH] . 
The electron dynamics in the presence of interactions in 
the dot [201 EI] and its spin/charge separation [22] have 




FIG. 1. Schematic view of the Anderson quantum RC circuit. 
Spin-degenerate electrons tunnel, through a quantum point 
contact, between a reservoir lead and a single level with a 
local interaction. An oscillating voltage Vq(t) is applied to 
the dot via a metallic gate of geometric capacitance Cg. 



been also studied. For small metallic islands the prob- 
lem has been addressed at intermediate temperatures [23] 
and in the many channel case |24j . In particular, the two- 
channel case has been argued to exhibit non- Fermi liquid 
behavior [25l [26] . Novel perspectives have been opened 
by recent experiments [27 30] where a significant dipole 
coupling between a microwave superconducting resonator 
and a quantum dot has been demonstrated. 

The low frequency admittance for the current / from 
the dot to the lead can be matched with the correspond- 
ing formula for a classical RC circuit 
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This allows to define a quantum capacitance Cq and 
a charge relaxation resistance Rq for the AC admit- 
tance of the system. This formula is related to the dy- 
namic charge susceptibility of the dot Xc(w) by the re- 
lation I{uj)/Vg{uj) — —iuje'^Xci^)- Xc{i — t') — i/h9{t — 
t') {[n{t),n{t')]) is the linear response function of the total 
occupancy n of the dot to a change in the gate voltage. 

These quantities have raised a large interest from the 
theoretical point of view starting with the seminal works 
of Biittiker, Pretre and Thomas [31311 [35]. In the quan- 
tum regime, the quantum capacitance Co provides infor- 
mation on the level structure |33j of the quantum dot. 
For a single channel in the QPC connecting the dot and 
the lead, Rq is universally fixed to /i/2e^ regardless of 
the QPC transparency. This prediction has been exper- 
imentally demonstrated [5]. It coincides with the lead- 
reservoir interface resistance j34j relevant in DC trans- 
port [351 [3S]. 

The universality of Rq still holds if interactions in the 
dot [37] or not too strong interactions in the lead [38] 
are taken into account in an exact manner. Increasing 
the size of the dot results in a mesoscopic crossover for 
Rq from /i/2e^ to h/e'^ [37 . For strong enough inter- 
actions in the lead, i.e. a Luttinger parameter below 
1/2, the system undergoes a Kosterlitz-Thouless phase 



transition to an incoherent regime where Rq is no longer 
quantized (38) . 

In this paper, we investigate the AC hnear regime of 
the quantum RC circuit where electrons carry a spin de- 
gree of freedom, as represented in Fig. [T] and the system 
is described by the Anderson model. Throughout the 
paper, we shall focus on the regime where the local in- 
teraction term U is much larger than the hybridization 
energy T such that charge fluctuations are small except 
at the charge degeneracy points, i.e. the Coulomb peaks. 
It includes in particular the Kondo regime where the spin 
on the dot is strongly correlated with the Fermi sea in 
the reservoir lead. Our analysis shall also include the 
more exotic SU(4) Kondo regimes relevant for dots with 
an additional orbital degree of freedom [39l [40] . 

The charge relaxation resistance of the Anderson 
model has been recently investigated by numerical renor- 
malization group (NRG) calculations [H], where it was 
shown that Rq develops a giant peak at zero temperature 
for Zeeman energies of the order of the Kondo temper- 
ature Tk- An analytical description of this peak has 
been given in the Kondo scaling limit |32]j on the basis 
of a Fermi liquid description valid at low temperature, 
in quantitative agreement with the NRG results. It also 
predicts the disappearance of the resistance peak at the 
particle-hole symmetric point, Ed = — C//2, where Sd de- 
notes the single-orbital energy on the dot. The Fermi 
liquid approach [15] is based on the identification of the 
low energy effective model, consistent with the Friedel 
sum rule. It allows one to derive a generalized Korringa- 
Shiba relation l44l 
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which relates the dynamical charge susceptibility Xci^) 
to the static ones Xa — ^d {n„) /ded. {n„) denotes the 
static occupancy of the dot for spin a. A similar relation 
was previously obtained for the spin susceptibility using 
the same Fermi liquid arguments [45] . Comparing Eq. (pi) 
with Eq. (fTl) , a general formula for the charge relaxation 
resistance 
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FIG. 2. a) Comparison of Rq as a function of the mag- 
netic field between NRG calculations (dots) (extracted from 
Ref. [IT]) and BA results (solid lines) for different Sd/U and 
U /T = 20, showing good agreement, b) Comparison of the 
universal function <^o{H/Tk) (solid) derived in the Kondo 
scaling limit in Ref. [42] with the function "I> (dots) derived 
by Bethe ansatz calculations from Eq. (Ill for U /T — 20 and 



close to the mixed-valence region {sd/U — —0.2) 
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is extracted. Xc = XT + Xi is the total charge suscepti- 
bility and Xm = XT ~ X^ the charge-magneto susceptibil- 
ity [32] . Xm should be clearly distinguished from the spin 
susceptibility, which is the derivative of the magnetiza- 
tion with respect to the magnetic field. The whole point 
of Eq. (l3|) is that Rq , a dynamical quantity, is expressed 
in terms of static quantities computable by Bethe ansatz 
(B A) . Deviations from universality occur in Eq. ([S]) when 
Xm 7^ 0, that is when both the particle-hole and the 
SU(2) spin symmetries are broken. 

In this work, we extend the Fermi liquid analysis to 
different parametric regimes by solving numerically the 
BA equations for the ground state [^5H15] and computing 



the static susceptibilities Xc and Xm appearing in Eq. ([3| . 
In the Kondo region, the robustness of the scaling form 
proposed in Ref. [32] is tested for finite parameters of the 
Anderson model, as shown in Figs.[2jb and[8] We confirm 
notably in Fig. [2Ja that the Fermi liquid result Eq. (Is]) 
agrees nicely with the NRG calculations of Ref. 41]. 
Out of the Kondo regime, departures from universality 
of the charge relaxation resistance as a function of the 
magnetic field were also shown within the Hartree-Fock 
approximation [2T]. Extending the BA calculations to 
the mixed-valence, empty orbital and valence-fluctuation 
regimes, we find that the peak in the charge relaxation 
resistance survives in these regimes, although its mag- 
nitude decreases in size with £d/U . Interestingly, even 




FIG. 3. Scaling form of the charge relaxation resistance Rq in 
the limits T <^ Sd <^ U (solid line) and F <^ U <^ Ed (dashed 
line). In circles, the resistance Rq is plotted for U — Ed- 



far in the valence- fluctuation region, i.e. for large £d/U 
and H c:i Sd, the peak is still present: Rq varies between 
/i/4e^ and /i/2e^ as a function of the magnetic field. The 
corresponding universal function for Rq, represented in 
Fig. [31 is derived analytically using perturbation theory 
and shown to agree with the BA calculations. In this 
region, the peak in Rq is not generated by breaking the 
Kondo singlet, but by the transition between different 
charge states of the dot. 

We finally give a further application of the Fermi liq- 
uid approach [HI |33] by considering an additional orbital 
degeneracy in the dot responsible for SU(4) Kondo be- 
havior at low energy 39, 40J. The existence of a Fermi 
liquid ground state [ISJ [SD] in the case of a SU(4) sym- 
metry allows us to derive an analog of Eq. ^. In the 
Kondo scaling limit, we predict, similarly to the SU(2) 
case, a giant peak in the charge relaxation resistance. 

The paper is organized as follows. Sec. [ll| explains 
how Co and Rq are calculated by solving the BA equa- 
tions for the ground state of the Anderson model once 



the Fermi liquid fixed point is determined. In Sec. Ill 
we study the range of validity of the Kondo scaling limit 
obtained in [32] ■ In Sec. IV we analyze the new scaling 



forms of Rq in the valence-fluctuation region. The peak 
in the charge relaxation resistance for a SU(4) symmetric 
Anderson model in the Kondo limit is presented in Sec. 
lYl 



II. FERMI LIQUID PICTURE 

The relevant model to describe the quantum RC circuit 
in Fig. [T] when the dot level spacing is sufficiently large 
and the transport is not spin-polarized, is the Anderson 



model [12132] 
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This Hamiltonian describes a single level, whose double 
occupation costs a charging energy U , weakly coupled to 
a non-interacting electron bath. The operators Cfco- and 
da annihilate electrons of spin a on the lead and on the 
dot respectively. The lead electrons are characterized by 
the single-particle dispersion relation Sk with a constant 
density of states vq. The total electron occupancy of the 
dot is n = n-|- -|- n^ with rio- — d^d^. The geometric 
capacitance Cg and the tunable electrostatic coupling Vq 
between the dot and the metallic top gate enter in Eq. (|4| 
through the interaction, or charging, energy U — e^ /Cg 
and the single-electron orbital energies Sda = —eVg — 
aH/2 where H is the external magnetic field, t is the 
amplitude for electron tunneling between the dot and the 
lead and we assume the hybridization constant F = tivqI^ 
to be independent of the magnetic field [5T] . 

It is a well established fact that the Anderson model 
behaves as a Fermi liquid at zero temperature [52, 53] 
for all values of the single-electron orbital energies £da- 
Moreover, the phase shift of quasi-particles at the Fermi 
energy is fixed by the dot occupancy through the Friedel 
sum rule 1541 
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For a time-dependent gate voltage £d{t) — e'd'^^uj cos(a;i), 
the form of the Hamiltonian follows from a quasi-static 
approximation [42l |43] , consistent with the Friedel sum 
rule Eq. ^ 
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where the a-^a operators describe quasiparticle states 
with a phase shift (5cr(£°) with respect to the original 
fermions Cka- 

The dot variables have disappeared from the effective 
Hamiltonian Eq. (l6]), although the memory of the dot 
is kept in the static charge susceptibilities Xo-- Practi- 
cally, the occupation number (n) = (n^) -I- (n^) and the 
magnetization (m) = ((n-|-) — {n^))/2 are static observ- 
ables and they are obtained by solving numerically the 
BA equations summarized in Appendix \K\ 

The study of (n) and (rn) identifies four regimes shown 
in Fig. |4^. A finite hybridization F between the dot and 
the lead smoothens the boundary lines between the differ- 
ent charge and spin states of the dot, as seen in Figs. |4|3 
and |4|:. The region where the charge is equal to 1 and 
the magnetization to 1/2 is called the local-moment re- 
gion. The transition to the empty (or doubly occupied) 
orbital regimes, where the charge is held fixed to zero 
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FIG. 4. a) Phase diagram of the isolated dot in the presence 
of a magnetic field, b) and c) are computed from the Bethe 
ansatz equations summarized in Appendix K] The domain 
Ed < —U/2 is deduced from the colored domain {ed < ^U/2) 
by particle-hole symmetry, b) Occupation number of the dot 
(n), for U/T = 20, reproducing the phase diagram in a). The 
boundaries between charge states are smooth functions. The 
smaller panel represents the static charge susceptibility Xc 
with its Coulomb peaks, c) Magnetization (m) of the dot, 
resembling the phase diagram in a), except at low energy. 
The difference is more visible in the smaller panel (logarithmic 
scale), where the local moment is screened below the Kondo 
temperature Tk (solid line). 



(or two), takes place in the valence-fluctuation region. 
The valence-fluctuation region is signaled by a Coulomb 



peak in the charge susceptibility Xc (visible in the smaller 
panel of Fig. l4]b) which defines the frontiers between the 
different Coulomb-blocked regions with zero, one or two 
charges. For Zeeman energies below Hi = rU/{U-\-2ed), 
the mixed-valence region is entered and the Coulomb 
peak deviates from the H = 2ed line touching the H ^ 
axis at e*^ = 0, where e*^ = Ed + F/tt ln(7re[//4r) is the 
renormalized orbital energy of the dot [IHIISI]. This devi- 
ation is presented in Fig. [5] The magnetization, shown in 
Fig. |4jc, shows a different behavior from the charge occu- 
pation of the dot. The transition line between a magne- 
tized and a non-magnetized state penetrates in the local- 
moment region following the Kondo temperature |47) 
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This is the signature of a strongly correlated ground state 
where the lead electrons screen the spin of the dot by 
forming a many-body Kondo singlet |55j . In general, 
this state cannot be described by standard perturbation 
techniques. In this paper, we circumvent this difficulty 
by solving the BA equations for the static quantities, 
combined with a Fermi liquid approach to access the low 
frequency behavior of the dynamical charge susceptibility 
Xc(w)- 



A. The quantum capacitance Co 

The quantum capacitance Co = e^Xc appears to lead- 
ing order in the frequency expansion of Eq. (IT]). The 
static charge susceptibility Xc can be calculated from the 
Bethe ansatz solution and is plotted in the inset of Fig. 
|4]b. It exhibits strong Coulomb peaks at charge degener- 
acy points for C/ 3> F, as a result of charge quantization. 
Xc is also represented in Fig. [s] as a function of the gate 
voltage for different values of the magnetic field. Fig. [5] 
illustrates in particular that Xc is insensitive to the mag- 
netic field until the Zeeman energy is of the order of F. 
In the Kondo region, the Kondo temperature Tr- is much 
smaller than F, and the peak in the charge relaxation 
resistance Rq thus develops in a region where the static 
charge susceptibility is independent of the magnetic field. 

When the Zeeman energy is above the hybridization 
constant F, the Coulomb peak starts moving following 
the H = 2ed transition line obtained for the isolated im- 
purity diagram in Fig. [4]a. In this regime, the Coulomb 
peak has a Lorentzian shape which can be derived an- 
alytically by just neglecting the spin down component. 
This procedure will be presented in Sec. |TVJ 

We stress that, in contrast with the non-interacting 
case, the quantum capacitance is not proportional to the 
local density of states as it is sensitive only to charge 
excitations and not to spin excitations. Hence, the Kondo 
peak in the density of states, which arises due to spin-flip 
processes, has no effect on the quantum capacitance Cq. 
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FIG. 5. Charge susceptibility Xc for U /T = 20. The circles 
(left to right) correspond to H/T — 0.1, 5 and 15 respectively 
and show the displacement of the Coulomb peak, also shown 
in Fig. |4]a. The solid hue is obtained for H/T = 0.0001 
and almost coincides with H/T = 0.1, showing the very weak 
dependence of Xc with the magnetic field in the Kondo regime. 
For higher magnetic fields, H/T = 5, 15, Xc converges to the 
Lorentzian form Eq. ( 19 1 (dashed lines) derived in the valence- 
fluctuation region. 



B. The charge relaxation resistance Rq 

The second term in the low-frequency expansion of 
Eq. Ilj describes the leading deviation from adiabatic- 
ity and introduces the response time scale RC to a slow 
drive of the gate voltage. Eq. (|2| derived in Ref. gS 03] 
gives the charge relaxation resistance Rq for all gates 
voltages and magnetic fields. Xc and Xm are both com- 
puted by solving the Bethe ansatz equation summarized 
in the appendix. Before discussing the results for Rq in 
the different regimes of parameters, let us note that the 
particle-hole symmetry of the Anderson model implies 
that Xm is an odd function of Ed + U/2 and thus van- 
ishes for Ed = —U/2 . As a result, the quantized value 
Rq — h/Ae^ is obtained at the particle-hole symmetric 
point irrespective of the magnetic field. 

In the Kondo region, Rq assumes the form 



^«-i^ 



^)Vo(.)a>o(^" 



(8) 



in the scaling limit [/ > E, e;^ > F, iJ < F. The function 
$o(a;) = xf'{x), plotted in FigJ2]b, is obtained from the 
universal form of the magnetization m = /(H/Tx) for 
the Kondo model [56 with the asymptotic behaviors: 




x< 1, 



a;> 1, 



(9) 



(Inx)^ ' 

where e is Euler's number. The function $o develops a 
peak when the magnetic field H is on the order of the 



Kondo temperature Tk- The envelope function 



(l-fy2) 
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depends on the asymmetry parameter y = l + 2ed/U and 
is shown in Fig. [8J It is obtained from the leading order 
charge susceptibility (insensitive to the magnetic field for 
iJ <C F) and from the derivative of the Kondo tempera- 
ture de^ InTR- = TTy/T. d^^ InT^ is an odd function of y 
such that Fo{y) vanishes at the particle- hole symmetric 
point y = in agreement with the above discussion. 

The robustness of the scaling form Eq. (IS]) for finite 
values of the different parameters of the Anderson model 
is discussed in the following Section. 



III. SCALING FORM OF THE CHARGE 
RELAXATION RESISTANCE 

By definition, the scaling form Eq. ([8]) is only an 
asymptotic behavior and it is of interest to evaluate how 
quantitative it is for real systems. In the general case, 
we extend the definitions of the two functions 
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such that they coincide with Fq and $0 in the scaling 
limit. In contrast to Fq and $0, F and $ do not depend 
solely on y and H/Tk but on all parameters of the An- 
derson model U, Ed, E and H . The range of practical 
validity of the scaling form Eq. (Isl) is tested below. 



A. The resistance peak in the function $ 

The departure of <& from $0 is studied in Fig. |6] by 
plotting $ as a function of the magnetic field H for dif- 
ferent values of C//F and the asymmetry parameter y. 
The Kondo temperature Tr- used to rescale the mag- 
netic field in Fig. [6] is obtained by numerically matching 
the low field behavior of $ with the expected asymptotic 
form ^q{H/Tk) ^ H/Tk for H <^Tk. The result for 
Tk is shown in Fig. [7] where it is compared to the Kondo 
temperature Eq. (It]) of the Anderson model. 

A first regime can be identified for C//F > 5 where the 
universal function ^q{H/Tk) is well reproduced in the 
Kondo region. The deviation between $ and $0 becomes 
sizable only close to the Coulomb peaks, where e*^ ~ 0, 
as seen in Fig. [5j At these charge degeneracy points, 
the peak in the charge relaxation resistance decreases in 
magnitude with Ed but does not disappear. The form 
of the resistance peak in the crossover from the empty 
orbital to the valence fluctuation region is discussed in 
SeclTVl 

For C//F < 5, Kondo physics is much less pronounced 
which results in a lowering of the peak in $. The agree- 
ment between the calculated Kondo temperature using 
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FIG. 6. Function $ Eq. (Ill for different Ed/U (squares and 
circles) and U /T calculated by the numerical solution of the 
Bethe ansatz compared to the universal scaling $o (solid line). 
The values of Tk for the numerical data are fixed by matching 
the values of $ at low fields to the linear behavior of $o Eq. 
Q. These are plotted in Fig. [t] 
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FIG. 7. Numerical values of Tk obtained (see main text) from 
the BA equations for U /V = 2, 5 and 20 (full circles, squares 
and empty circles respectively). They are compared to the 
analytical formula Eq. ([7| (solid lines). 



B. The envelope function F in the Kendo region 

Fig. [5] demonstrates that, as long as one remains in the 
Kondo region, the dependence of Xc on the magnetic field 
can be safely neglected. In Fig. |8]a, our BA calculation 
for the F function at zero magnetic field, represented by 
the dashed line, is in very good agreement with the NRG 
data extracted from Ref. 01]. It remains however far 
from the asymptotic function F(j{y) even though U /T = 
20 in Fig.jSJa, see also Ref. gH. 

The convergence of F to the asymptotic form Fo{y) 
as a function of U/T is illustrated in Fig. Mb where it is 
shown to be slow. A more quantitative analytical expres- 
sion for F can be derived by including the next to leading 
order corrections to the charge susceptibility, namely |43| 
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in Eq. (111. The result for F is the function i^o(y) with 
additional T /U corrections. It is in much better agree- 
ment with the BA calculations and the NRG results from 
Ref. [41] than _Fo alone, as shown in Fig. [SJa. 



IV. THE VALENCE-FLUCTUATION REGION 



our fitting procedure and Eq. (l7| is also degraded as 
shown in Fig. [71 



The meaning of Eq. ([Sl is restricted to the Kondo region 
where a Kondo temperature can be defined. 



As we already saw in Fig. [6J the peak in the charge 
relaxation resistance decreases in magnitude at the edge 
of the Kondo region, in the mixed- valence region around 
ej ~ 0. Below we discuss the fate of the resistance peak 
as Ed is further increased to explore the empty orbital re- 
gion Ed S> r, and the valence-fluctuation region at higher 
magnetic field. As we shall see below, the resistance peak 
does not disappear although its magnitude does not scale 
with U/T in this region. 

The peak in the charge relaxation resistance can be 
derived analytically in the regime e^ 3> F by standard 
perturbation theory. In this regime and for arbitrary 
magnetic field, the two states of the isolated dot forming 
the low energy sector are |n = 0) and |n = l,t) as shown 
in Fig. [9J The absence at low energy of the spin down 
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FIG. 9. Spectrum of the dot isolated from the lead on the 
left. For a positive Sd in the presence of a magnetic field, only 
the states |0) and |l,t) compete in the low energy sector. 
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FIG. 8. a) Comparison of NRG (points extracted from the 
results of Ref. i41j), BA (dashed line) and analytical (solid 
line) results for the function F with U/T — 20. b) Approach 
to the scaling limit Fq Eq. ([To| for different U/T. The dotted 
lines are obtained by BA while the solid ones correspond to 
the perturbative result Eq. (12 I. 



component implies that the ground state does not exhibit 
strong correlation and can be described analytically using 
perturbation theory. The unperturbed Hamiltonian is 
obtained by setting the tunneling involving spin down 
electrons 



t^ {(^ildi + dl'^H 



(13) 



to zero. In that case, the number of spin down electrons 
on the dot is a constant of motion and the Hamiltonian 
can be diagonalized separately for n^ = (low energy) 
and n^ = 1 (high energy). It gives an exactly solvable 
resonant level model 



H' = ^£fc4.^Cfct+(edt 
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(14) 
for which the charge relaxation resistance is h/2e^. 

Let us call \ipo) the unperturbed ground state, with 
n^ — 0, characterized by the spin up electron occupancy 
on the dot 
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The perturbation due to the tunneling term Eq. ( 13 1 
gives the first order correction to the wave function 



i^i)=^e( 



1 



Ek-U -Ed- 



j^dlckin^+ 



(16) 



Efe — £d 



-jfdlckiil - n^) \tpo) 



The projectors n-^ and (1 — n^) are necessary to determine 
the part of |V'o) with a spin up electron on the dot and 
the part with no electron. This implies the presence or 
not of the interaction energy U in the denominator of 
Eq. (16 1. The values of the spin a populations for the 



corrected ground state |i/)o) + \tpi) are 
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corresponding to the static susceptibilities 
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X; 



1 - {n^f 



(%> 



(18) 



(£d + f + C/)2 



X? 



ed + f + (7 




We have introduced 



X? 



r 1 



(19) 



the spin up susceptibiUty in the absence of the spin down 
component. 



The static susceptibihties of Eq. ( 18 ) are combined to 
give Xc = Xt + Xl and Xm = Xt - Xl- Substituted in 
Eq. ([3]) , they give an analytical expression for the charge 
relaxation resistance R„ which still exhibits a peak as 



a function of the magnetic field, as shown in Fig. 10 



Fig. 10 also compares the analytical expression for Rq 
with the BA calculations and shows an excellent agree- 
ment already for e^/r = 6. The peak height occurs 
around h/2e'^ and for H ~ 2ed- At this point, the spin up 
charge fluctuations are maximum, see Eq. (19 1, because 
the states \n = 0) and \n = l,t) are degenerate for the 
isolated dot when H = 2ed, and the spin down fluctua- 
tions remain small. Hence, the resistance is around h/2e'^ 
as in the single- channel spinless case. The position of the 
maximum of the resistance can be found perturbatively 
from the analytical solution Eq. (Il8|) 



H ^^ F UjAed + U) 
2ed TTEd (2£rf + t/)2 ' 

— (l ^ ^ 
'2e2 1^ ^ 7red2ed + U 



Ra 



(20) 
(21) 



The analytical expression obtained for Rq from 



Eq. ( 18 1 can be further simplified in the limit £d ^ F. 
For X — H/2ed < 1, the universal form 



Rn 



h 



1 



4x"^ 



(a;2 



1 



(22) 



is obtained. This result is independent of V because the 
unperturbed ground state is |0) when F is sent to zero. 



'? 4? 




HjlZd 



FIG. 10. Comparison between Kq obtained from the analyti- 



cal results Eqs. ( 18 1 (solid line) and the numerical solution of 
the BA equations (circles) for U jV — 20 and Stj/F — 6. 



The doubly occupied state is therefore reached only to 
second order in perturbation theory and can be neglected 
to leading order. For x > 1, the unperturbed ground 
state is 1 1, t) for a vanishing F and the form of Rq depends 
on the ratio Ed/U. For Sd ^ U, we recover essentially a 
non-interacting resonant level model and the resistance 
is also given by the universal form Eq. (22 1 for x > 1. 



For Ed ^U however, the charge relaxation resistance Rq 
is frozen to h/2e^ for all x > 1. Both these universal 
limits are shown in Fig. l3] The reason is that the doubly 
occupied state is forbidden for infinite U such that the 
spin down states cannot be reached within first order 
perturbation theory. Hence we only have spin up charge 
fiuctuations, xi ^ 0, and we recover the universal result 
of the spinless case Rq = h/2e'^. 



V. THE SU(4) KONDO CASE 



We extend our discussion to the more exotic case of 
a SU(4) Kondo effect [53] ■ This situation is relevant for 
certain quantum dots with an additional orbital degree of 
freedom that is conserved during lead-dot tunneling pro- 
cesses [57] . For example, ultra-clean carbon nanotubes 
have a natural orbital degeneracy that arises from the 
clockwise and anti-clockwise motions of electrons around 
the tube. We label here the orbital index by / = -I-, — . 
The model has now four transport channels in correspon- 
dence with the four available single-electron states in the 
dot: |-f,t), \+,i), |-,t) and |-,i). We label these 
four states by a quantum number r — 1, ... ,4 respec- 
tively and use the same index for the conduction elec- 
trons in the lead. The Hamiltonian takes the form of a 
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SU(4) Anderson model 



H ^ 



kr 



kr 
t<t' 



-dl 



rCkr 



(23) 



where the meaning of the operators and notations are the 
same as in Eq. El). For temperatures much below the in- 
teraction energy U and T <^ U, the charge on the dot is 
frozen to 1, 2 or 3 depending on the gate voltage Ed- Per- 
forming a Schrieffer- Wolff transformation 58J , one finds 



HSW = X! Sfc^rCfer + Wq^ Ar^k'T 



kk' 



2 '^ 

kk'TT' 



4rCfe' 



dl,dr- Isj^TT' 



(24) 



where q denotes the dot occupancy in the low energy 
sector and A^ == 4. The generalization to any N and q 
is straightforward. The values of the potential scattering 
and the Kondo coupling constants Wg and Jg are given 
by 



1 



1 



A/ 


\ed + {q-l)U Ed + qU 
e / q N-q " 


^q 


N [ed + iq-l)U ' Sd + qU^ 



(25) 
(26) 



The potential scattering term vanishes for Sdwo = (1 ^ 
q — q/N)U. An exact mapping to the SU(N) Kondo 
model [59] is then obtained 



HsU{N) - 22 ^''^kT 



^kr 



j'gS ■ T, 



(27) 



where Jg 



2t^ 



N^ 



jj- g(pf_g) ■ We switched to the basis of gen- 
erators of SU(N) [55ti61j . such that an anti-ferromagnetic 
coupling between the spin S = ^^ j+ ^ 

Z-tkk'TT' ^kr ^ 



— dr' of the im- 



V d^ - 

purity and T = ^Ylikk'Tr' '^kr^^^'^k'T' of the lead is made 
explicit. A is the vector composed of the A^^ — 1 matrices 
which compose the A x A fundamental representation of 
the SU(N) group. Their explicit expression in the SU(4) 
case can be found in [62^ . 

As mentioned in the Introduction, the low energy 



fixed point of the Hamiltonian Eq. ( 24 1 is a Fermi liq- 



uid and the Fermi liquid approach 42, 43 introduced 
in Sec. [ll] is also applicable to this model. Defining 
Xt — —d (tIt) /dsd as the r-dependent static suscepti- 
bilities, the charge relaxation resistance is found to be 



Ra 



2-^r Xt 



2eME.Xr) 



(28) 



The emergence of logarithmic singularities prevents the 
study of the Xt susceptibilities by perturbative methods 
below the SU(4) Kondo temperature jH] 



T'* = X'e'^^'-^""'^'''' 



(29) 



where V ~U^£d is the effective high-energy cut-off of the 
model whose precise form is not needed here. 

Following the line of reasoning developed in Ref. [42] , 
one can derive the behavior of Rg in the presence of a 
magnetic field. We first switch to a more convenient basis 
that separates the charge, spin and orbital degrees of 
freedom, namely 



(30) 



In addition to the total charge susceptibility Xcj we have 
introduced the charge magneto-susceptibility Xm, as in 
the SU(2) case, and its orbital counterpart, x^, which 
measures the sensitivity of the orbital magnetization to 
a change in gate voltage. Xmv is obtained from the differ- 
ence between the spin magnetizations of the two orbital 
states. 

Substituting the new susceptibilities in Eq. ( 28 1 , the 
charge relaxation resistance is found to be 



Xc \ 


/I 


1 


1 


1 \ 


(XI 


Xm 


1 


-1 


1 


-1 


X2 


Xv 


1 


1 


-1 


-1 


X3 


Xmv / 


\1 


-1 


-1 


1 / 


\X4 



Ra 



Xr. 



X.V 



AVT. 



(31) 



the analog of Eq. (|3| in the SU(4) case. At zero magnetic 
field, the spin and orbital degeneracies are not broken 
such that Xm = Xv = Xvm = and a universal resistance 
Rg — h/Se^ is obtained. At finite magnetic field, only 
the spin degeneracy is broken and Xv = Xvm = 0. 

In the limit U ^ T and for magnetic fields of the order 
of the Kondo temperature Eq. ( 29 ) , the magnetic field de- 



pendence of the charge susceptibility Xc can be neglected. 
Assuming that the results of Cragg and Llyod [B3] are 
also valid in the SU(4) case, such that the leading po- 
tential scattering term in Eq. ( 24 ) is unaltered along the 



Kondo crossover, the Friedel sum rule leads to 



Xc ^4:Vode^W„ 



q 4-g 

iSd + iq-l)Ur (Sd + qU)^ 

(32) 
in the sector with q charges on the dot. 

As in the SU(2) case, the form of the charge magneto- 
susceptibility Xm can be derived from scaling arguments. 
In the Kondo limit, the magnetization of the dot m — 
X)/ CT '^ {"^iry) /2 is a smooth and monotonous universal 
function fg(H/Tj^) that starts at zero at vanishing mag- 
netic field and saturates at 1/2 (resp. 1) for large mag- 
netic fields, when g = 1,3 (resp. q = 2). Differentiating 
the magnetization with respect to Ed, one obtains 



Xra^2d,JnT^^<Pg{H/T^^) 



(33) 



where we defined the universal functions ^g{x) — xf'{x). 
From the general form of the functions fg{x), we expect 
that the functions ^q{x) have a similar peaked shape as 
the function $(a;) of the SU(2) case. Using the expression 
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FIG. 11. Envelope functions Fq in the sectors where the 
charge is frozen to 1, 2 or 3, that is when Sd/U £ [-1,0], [- 
2,-1] or [-3,-2] respectively. The function becomes zero in the 
middle of the Coulomb valleys, while the circles correspond to 
the values of Ed/U for which the potential scattering couplings 



Wq in Eq. ( 26 1 become zero 



Eq. (29) of the Kondo temperature, we obtain to leading 
order in T /U 



IT 2ed + {2q-l)U 



2r 



u 



$„ 



H 



(34) 



where the prefactor 2ed + (2g — 1)U essentially comes 
from the derivative of the Kondo temperature Eq. ( 29 1 . 
Combining Eqs. (32) and (34) into Eq. (31), we find a 



scaling law in the Kondo limit 



7? - ^ 



U 



^g(y,)* 



H 



(35) 



similar to the SU(2) case. Thus a giant peak in the 
charge relaxation resistance, proportional to (fJ/T)^, also 
emerges for a SU(4) symmetry. The envelope functions 



Fq{y,) 



{vl - 1) 



[1 



vl 



vM-'^W 



(36) 



depend on the charge q and on the variable yq = 2edlU + 
2q— \. 2/g is defined such that j/g = ±1 at the Coulomb 
peaks and y^ = in the middles of the Coulomb valleys. 
The envelope functions corresponding to the three charge 
sectors q — 1,2,3 are represented on the same plot Fig. 
11 as a function of £d/U . 



Interestingly, the function F2 coincides with the SU(2) 
function F up to the multiplicative factor 16. Instead, 
in the sectors q ^ 1 and 3, the envelope function is 
asymmetric, which gives an experimental signature dis- 
tinguishing SU(2) and SU(4) symmetries. We also notice 
that the values of £d/U , for which the envelope functions 
Fi^ vanish, do not coincide with the locations of zero 
potential scattering, i. e. Wq = in Eq. ( 26 ) , represented 



by circles in Fig. 11 



We expect that the approach to the 
Kondo scaling behavior is faster at those latter points 
since they are free of potential scattering and exhibit only 
Kondo coupling. In addition, the envelope is close to its 
maximum at these points, in contrast with q = 2 and the 
SU(2) case where the envelope vanishes as imposed by 
particle-hole symmetry. 



VI. CONCLUSIONS 

In this paper, we performed a thorough study of the 
quantum capacitance and the charge relaxation resis- 
tance for the Anderson model. We applied a Fermi liquid 
approach, where the low energy effective model is derived 
consistently with the Friedel sum rule, that allowed us to 
express the charge relaxation resistance in terms of static 
susceptibilities. The susceptibilities are computed from 
the Bethe ansatz equations describing the ground state of 
the Anderson model. The accuracy of our approach was 
tested by comparing our results to NRG calculations [^ 
or perturbative calculations both in the Kondo and in the 
strongly asymmetric regimes. The analytical predictions 
given in Ref. [H] for the peak in the charge relaxation 
resistance are shown to apply in the whole Kondo region 
for C//r > 5. The persistence of this peak was demon- 
strated in the valence-fluctuation region, both numeri- 
cally and from a direct perturbative calculation. More- 
over, we showed how the Fermi liquid approach can be 
extended to the SU(4) symmetric case where a similar 
peak emerges in the charge relaxation resistance. 

Overall, this work constitutes a specific and detailed 
example of how the effective Fermi liquid theory can 
be used to derive the low frequency dynamics of quan- 
tum impurity systems. This does not include, of 
course, systems and regimes in which non- Fermi liquid 
physics [25l [26] dominates such as impurity models with 
overscreening. Further extensions of this work could in- 
clude the study of non-zero temperatures and higher fre- 
quencies [411 I64j where inelastic processes play an in- 
creasing role. 

We acknowledge T. Kontos for useful discussions and 
thank the authors of Ref. IT for providing us their NRG 
data. KLH acknowledges support from DOE under the 
grant DE-FG02-08ER46541. 



p(fc) + g\k) I dpp{p)R[g{k) - g{p)] + g'{k) / dX<j{X)s[g{k) - A] = 5''(fc), 
rQ 



a{X) 



dA'cr(A')i?(A - A') + / dkp{k)s[X - g{k)] = S^iX), 
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Appendix A: Bethe ansatz equations for the ground state of the Anderson model 

A striking feature of one dimensional quantum systems [65^ is the possibility to have a separation between charge 
and spin degrees of freedom for electrons at low temperature. In the case of the Anderson model, spin and charge are 
carried by different excitations called spinous and holons respectively. Their densities of states are denoted p and a. 
They satisfy the following Bethe ansatz integral equations [JBHiS] (we follow the notations of Ref. [iS] ) : 



(Al) 
(A2) 

(A3) 
(A4) 

(A5) 
(A6) 



with the source terms given by 



1 



SP{k) = — ^l + g'{k)J dpR[g{k)-g{p)] 



L 



A{k)+g\k) dpA{p)R[g{k) - g{p)] 



/CX3 
dks{X-g{k)) 
-CX3 



1 , A(fc) 
2tt L 



We have introduced the functions 

1 



Rix) = 



2n 



duj 



1 + el'^l' 



m 



k-eg- U/2 
2UT 



s{x) 
A{k) 



1 



2cosh(7ra::) ' 

r 1 



^(fc-e<i)2+r2- 



L is the size of the system and the holon and spinon densities can be split in a conduction and impurity (dot) part 



p{k) = p,{k) + 



Pijk) 
L '■ 



ct(A)=ct,(A) + 



a^{X) 



(A7) 



The linearity of Eqs. ( Al I and ( A2) implies that the conduction and impurity terms decouple. The former fixes the 



macroscopic properties of the system, i.e. the global magnetic field H and the position of the valence level Ed, 

(•B 



H 

2^ 



dkpc{k), 



1 / U 

^h+2 



dAcrc(A), 



while the latter gives the occupancy (n) and the magnetization {ni) of the dot, namely 



dkp.i{k), 



{n) = l 



dX(Ti{X). 



(A8) 



(A9) 



These equations hold exclusively for Ed > —U/2 and H > 0, while the results for Sd < —U/2 are obtained by 
particle-hole symmetry. 

The zero magnetic field case H = and the particle- hole symmetric point Sd = —U/2 are obtained by setting B 
and Q respectively to — cx). In these cases, the BA equations for p and a decouple and an analytical solution can be 
constructed on the basis of the Wiener-Hopf method [47 . 
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